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Abstract 

A general method for eliminating the angle variable from the 
equations of a perturbed periodic motion and for deriving an "adiabatic 
invariant" J has been given by Kruskal, and for a special clans of 

t 

2 

Hamiltonian systems, ilclJamara and Whiteman have shown ^ to order ) 

that J is related to a set of invariants I obtained from the expan- 
sion of Poisson bracket relations. In this work, an order-by-orde~ 
algorithm for Kruskal's method is introduced and a new 3et of invariants 
Z^ is obtained. It is shown that these invariants bear a close relation 
to those obtained from the Poisson bracket expansion, and in the special 

case investigated by Mchnm&ra and Whiteman, the relution between I and Z. 

t x 

may be brought to the same form as the relation between 1 and J 
derived by those authors. Finally, the relationship between and J 

is examined and arguments are presented that in certain cuse3 the two 
i . re equal to all orders. 
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IilTkODUCTlOi* 

Let a perturbed periodic mechanical system be given , described by n 
canonical variables collectively represented by the vector y and by a 
Hamiltonian H dependent on a small parameter , of the form 

H = y 1 + £ £ k H (k, (y) (l) 

• 

(here and in what follows, superscripts in parentheses denote order in £- ). 

The vector y may be viewed as the solution of the Hamilton-Jacobi equation 
for the unperturbed motion, yielding an action variable y^ , its conjugate 
angle variable y & and a set of : ther variables y^ which are constants of 
the unperturbed motion. 

A perturbation expansion may now be employed to eliminate y^ and y^ 
from the equations of motion to any desired order in & . One such technique 
has been devised by Kruskal^ and will be described in more detail further on 
(this method is also applicable to non-canonical systems). By Kruskal's approach, 
a near-identity transformation from the variables y to new "nice" variables z 
is performed, so that of the n first-order differential equations describing 
the evolution of z , (n-l) may be separated and solved independently as an 
autonomous set. furthermore, using these "nice variables", an "adiabatic invariant" 
J , which is a constant of the perturbed motion, may be expressed to any order of E • 
If J is used to eliiuinate , one winds up as required with a mechanical system 
containing only (n-2) independent variables. 

An alternate method of deriving an invariant I for the system described 
by eq. (l) is based on an expansion of the Poisson bracket relation 


L i. h] 


0 


- 2 - 


(*) 


l 2) 

This expansion has been described by Whittaker and was further explored 
by i-'cIJanara and »vhi toman ^ ^ henceforth referred to as HchW). The latter 

2 

authors were able to show — by a lengthy calculation anu only to order t. — 
that the invariant thus obtained is, in a special oaso, related to kruskal's J . 


In what follows we shall show that for a system given by eq. (l), nruakal'3 
method may be modified to yield a family of invariants in a completely different 
way than that used in deriving J . It will then be shown that ~oese invariants 
are solutions of eq. (2) , are related to the invariant of KcNW and also are 
connected with the adiabatic invariant J • 


DOTATION 

In order to obtain concise expressions, the notation used here departs 
somewhat from that of Kruskal. In kruskal's work, y and z_ have (n-1 j compo- 
nents and special symbols V and ^ stand for what we here uenote by y R 

and z . In what folLows, such vectors with (n-l) components, excluding the 
n 

angle variable, will be denoted by a tilde, e.g. y , zf . Furthermore, if the 
canonical variables in y are arranged in order 

y = ( p . q ) (3) 

we shall define a conjugate vector 

7 

so that y x = y n . y R = -y x , 


( q » -P ) 


U) 


The use of y enables one to write Poisson brackets concisely (the sumiuation 
convention is henceforth used in all sumuatipns over n or (n-l) components) 
as 

La, b} - (^a/^y.)('Db/Dy i ) (b) 

and Hamilton's equations become 


dyj/ dt (6) 

Finally, 

Y've shall assume (as in Kruskal's work) that the basic period in the dependence 
on y n is unity and denote quantities averaged over y R (which clearly depend 
on y only) by angular brackets 

a > f a dy 

J 0 n 

KRUSKAL'S BXPAK3IQK 


Kruskal's method does not require the system to be canonical but assumes the 
evolution of y to be given by equations of the form 


#Q 

dy / dt = I! g k ^(y) (7) 

k=0 

(k) 

where the components of g are periodic in y^ with period unity and where 
g^ has only one non- zero component, namely the n-th. Since wo are interested 
in systems for which (/) reduces to (6) with H given by (l), we shall assume 
that this component is unity 



( 0 , 0 ••• 0 , 1 ) 


(B) 


although what follows can be extended to more general cases. We now soek a 
near-identity transformation to new variables 
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♦ Z t. k $ U) (y) 19) 

k=l 
satisfies 

r k . (k) ,*v\ , 0 

A t h ^z) UO) 

k=0 

does not appear on the ri t ;ht-hanu side, the 

first (n-l) equutiono of (10) form an autonomous system which may be solved 

independently, the solution then being substituted in the remaining equation 

to provide the evolution of z . 

n 

Substitution of (y) in the l.h.s. of (10) gives 


dz, / dt = dy / dt 

♦ L e. 

fcsl 

k OJ^Zy^Uy/ dt) (li) 

- Zfc k 

k=0 

t 

* i»| lk - J, /*y t ) } 

- L f. k 

k=0 

, iu 

- Ci^ (k) Ay n ) - Z2 g U, -v$ ik ~ 

o-i 

with the V operator defined in 

y space. 

.of (10), 

We now convert the r.h.sUo depend 


on y as well, using the Taylor expansion operators derived as follows, 
refining the exponential of a differential operator by means of the power 
series for e , one can formally express the Taylor expansion of h v '(][) 
as follows ( * stands for "operates on") 

h tk) (|) = h (k) ( y + 

*• ei? ( . h U) (y) 

= Z t 3 tS k \y) U2) 

j=0 " ” 


such that the evolution of z 


dz. / dt 


with h^= g^ • Since z 
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where the are differential operators obtained by exp/inding the expo- 

(4) 

nential. The first few of thorn are ' 


,(0) = 

1 



d) 

v 






(13) 

12) = 

S (2, -v 

+ 7 ^5 (l) t77 


(3) 

% Vi) - v 

- S (1 ) t u) «vv 

♦ (l/ 6 ) V’V 


Because T^ equals unity, one can separate the h ^ term from the rest. 
The r.h.s. of (10; then becomes 

J_ f.“ |h (k) (y) + "z T (J) . h U ’ J, (y) j (14) 

Both sides of (10) are now functions of y and therefore the equality holds 
independently for every order. One thus obtains a set of equations 

"SS (k V'Sy n - h !k) (y) = * (k) (y) (15) 

where 

^ kJ (y) = £ T (J) . h (k - j) - g (l ° - £ g (j) . v< (k ' j) (16) 

- ' ' J=1 " 

and 

l 

k — 1,2, ... 

This may be used as a recursion relation. Suppose all quantities entering 

(k) i 

here are (like the g v ) either periodic in y^ or independent of it. Then 

(kJ (k) 

will possess a "secular” part < x ) indepenuent of y^ and a 

purely periodic part averaging aero 

( A 1 * 0 )^ = > (k) - 
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On the l.h.s. of (lb), is wholly secular by definition, while the 

1 ^ ) 

other tern there is purely periodic, since any secular part of £ is 

(k) 

removed by the differentiation. Thus if ' in given w» can derive 
the k-th order quantities through 


h tk \y) = - ( (17) 

^ U) /T>y n = - < ,V kJ > lie) 

from which 

£ k) ■ I 0 n {^ (k) ■ <^ (k) > j % * U9) 

with an arbitrary wwoular vector. In Krusknl's work, /i^ k is 

chosen to vanish, so that at y = 0 the vector z is identical with y ; 

n — 

we shall denote the functions thus obtained by . If however we only 

Ik) 

demand that z, should satisfy equations of the form (10), J* may be 
arbitrarily chosen. In what follows we shall make use of this free choice 
in order to endow the "nice variables" z with additional desired properties. 


THE NEW INVARIANT 


So far we have treated the general case of kruskal's expansion (apart 
from our choice of g^°M with no reference to the canonical character of 
y . Taking this now into account, one finds from (1) and (6) 


g^ k) * - ( 20 ) 

Defining for convenience 



and using (16) and (b) ; gives 
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(k) 


n 1 T (JJ * h U “ J \y)l + £ [h^.^ J Ul) 

j=»i L ' j j-i 


We now pose the following question; la It possible, by proper choice of the 
In (l9)f to make some component z^ of £ a constant of the 

motion? 

(k) • 

If 13 conserved, this means that h^ vanishes for all k , which 

Ik) 

in turn implies the vaniB>.ing of ( ' y . tty the last equation, this 

reduces to 


§ < [■“• sM > 


( 22 ) 


U) 


Suppose that at the stage when eq. (22) is reached the ‘ have been 

derived up to and including the (k-2) order. Ve can then fulfil {22) by 
(k-l) 

choosing yi-' to satisfy 

< k U) >] = {[**>. %?>]) * t {[^K $«-»]) 

1 j “2 

where the r.h.s. is as:ju’ued to be known at that rtage. The above equation is a 

linear first-order partial differential equation, and solutions in general do 

exist. Deriving them explicitely is another problem, however. For the special 

case when all ti v 1 with k >1 vanish, McNamara ar*d Whiteman (who arrive at 
obtained from first principles . a llow, / . , 
a similar equation) V formulas which to be derived (fcr i=l , 

which is the relevant cose, as will be seen) up to k=3. 

A different approach to the problem will be outlined in 

the last section of this article. 

The iteration for can thus be carried on — provided it can be 
started. The first additive function encountered is » U8 ®^ in ensuring 

the vanishing of { • Thore exists no adjustable variable to ensure the 


(23) 


- e - 


the vanishing of <i ) , so the iteration can be started if and only 

i 

if this tens vanishes of its own accord, which in turn implies 

< , W 1 V'-|>y i > - 0 (24) 

If (as assumed) y^ ent> s U only through angular tenon, this will 
certainly hold for i-1, since the y derivative vmich is applied in that 
case removes the secular part of H^ / , leaving a purely periodic function. 

Thus an invariant Z^ of the type discussed here may in general be derived. 

( 1 1 

If H v ' itself is purely periodic, other invariants may be 
generated for i ^ 1 . 


THK PQIJJQ1J dKAClJoT METHOD 

(3) ( 2 ) 

McNamara and Whiteman , following Whittaker , derive an invariant I 
(in their notation: J) in the following way. Lot I have an expansion in £ 

1 = £, £. k i (k) (y) 12!>) 

fe=0 

and let H bo expanded as in (l) (this is a slight generalization: in the 
cited work, all li ^ with i>l vanish^;. Then 

t 1 . «] = 0 = Z] if' jC 126) 

Hy (t>) 

[l (kj , h(°)] = [l u) , y x l = ^I U V^y n (27) 

So the equation for the k-th order is 

k-1 

'dl^V^y » - Hi (2b) 

n j*l 

(y) v 

which defines I * recursively within an arbitrary function of y 
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1^ = - f JZ L l(j) » dy„ ♦ G^ k '(y) (29) 

J 0 j=° 

- l< k > ♦ G U) 

> (k) 

By (l’j) and (21), this is exactly the same as the equation for ' , 
vk) 

provided all h^ vanish. 

Consider again eq. (20) : due to the derivative* its l.h.s. will 

be purely periodic, but unless special steps are taken the r.h.3. may well 

contain a secular part, tfe therefore must assume that at the stage at which 
(k) (k-l) 

I is being derived, G has not yet been determined. The r.h.s. can 

then be mad”' purely periodic bj requiring 

[c (k " 1; ,<H (l) >] = H ]) " (30) 

which has the same fora as (23). 

The lowest order of ^26) gives 

[l<°\ u (0) ] =0 = ^ °>ny n 

so that 

l(°) . X (0)(J) 

The iteration can be started only if 1^^ makes the r.h.s of (28) purely 
periodic for k=l : 

[V 0j , < = 0 (31) 

Obviously, any function of can be chosen as 1^ ; McNV choose 


l(°) = < H (U> 


(32) 
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This however ia not entirely satisfactory, since we 

<0) 

exvect I to tend to some natural invariant of the unperturbed system in the 
limit of vanishing £. , independent of the choice of . A more suitable 

choice is 

T (0) 


- yi 


03) 


which also aatial'ies (3l)» since (24 ) holds for i-1 . With this choice, I 

equals the action variable y^ in the unperturbed limit, a property shared by 
the invariant previously derived and also (it may be shown; by hruskal's J 

The alternative choice, made by HcU'V, will be explored in the next section. 

With 1^ chosen as in (33)* it also equals and it is easy 

to 3how that the expansion equations of I match those of the invariant Z. 


stage by stage. On may then match 


•00 




(k) 


0 (k) _ ^(k) 


( 1 ) 


by making identical choices of the arbitrary functions of H ^and by virtue 


of (33) satisfying (3l)» of y^ as well) which can be added to G' and 


to 


at every stage. 


A 

THE CHOICE 1^ ~ 


• ( 0 ) 


fccMaiaara and Whiteman chose I v ; as in (32), in the special case where 
(k) 

H Vanishes for k > 1 . In that case, their recursion continues with 

y. 


lU) ■ - J 0 “ [<* u) >, > u) . 3 K 


. 0 ) 


3(1) 


,( 1 ) 


(34) 
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jV 1 ). <H (1 VJ = <[h (i) , 1 u, ]> Ob) 

Sxpanding the invariant by Kruskal's method, lor the same H, 

^^1) . ^H (l) >Oy n 

and by (19) 

- H^\y n ^O) 

Fron (23) and the preceding equation, noting that secular iunction3 may be 
taken out of the averaging brackets 

jV^. <« ll) >] = - [<U (1J >. H (l, (y n =0)] 

Let us select 


1 (1) 
A 

= - <(h^ 1j> > + 

H U )(y n *0) 

so that 



S 1 ’ 

I 

H 

II 

|h)> 

Thie second order equations then are » by (2l) 

( .) 

3 1 

= - /“ [“ U) * 
J o 



[ r ;«.<•">] ■ <[•“’. i; M l> 

This demonstrates that, with matching choices of aduitive functions 

x w = .^ 2 > 

1 

c h) 
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If additive functions of higher orders are matched as well, the subsequent 
parts' of and of I are identical, except lor an extra order of £. 
in the former. That means 


- EA I - l^° ; ) = z - ^ 

1 1 1 


or 


H ■ n 


( 56 ) 


2 

a relation resembling that obtained (to order ) by McftW , except that 

in their result J appears in place of . 


RELATION TO J 


Kruskal defined the adiabatic invariant J by 

1 


- fc 


c k d \ 


j 2- Pk dz n 


( 37 ) 


where the integration is carried over a set of points ("ring") sharing the 

same £ and differing only in z . For details about J , the reader is 

referred to Kruskal' a article^; its value is independent of the canonical 

set used in its derivation, though the components of y are the best choice 

for this role, if they form a canonical set (the inverse transformation z -* y 

must then be derived by means of the txpinsion operators 113) )• Here we shall 

merely sketch out the connection between Z, and J without deriving the 

t 

details. 

Suppose that among th'i many sets of "nice" variables possible, differing 
« (k) 

in their choices of j * but all obeying (10), there exists a set (or a 
family of sets) that is canoidcal, with z^ conjugate to z R . Obviously, 
this set, too, can be used in deriving (37)* leading immediately to 
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The existence of nice canonical variables has been proved by kruskal^ in 

(k) 

appendix 2 of his article. It is furthermore possible to express the which 

generate such variables. The details of this derivation are somewhat involved and 

will therefore be described in a separate article; here we shall just assume that 

(k) 

these MV' are known. Then the transformation which they describe belongs to 
the (much larger) frmily of transformations which make a constant of the 

motion, each of which in its turn provides a solution to the corresponding 
Poisson bracket expansion. Thus among the many possible solutions of and 

I, there exist such ones for which 


I = z 1 = J 
(k) 

Practically, given the M : which make z_ canonical, these functions 

offer ' 

probably'fthe best way of deriving or I , since they are known to solve eq. (30), 
The only other possibility for solving this equation is to use the formulas of 
McNW, which are valid for the lowest few orders only and are specifically 
related to the choice of 1^ given in (32). 
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